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Abstract 

Platinum is widely used as a pressure calibration standard. However, the established thermal 
EOS has uncertainties, especially in the high P-T range. We use density functional theory to 
calculate the thermal equation of state of platinum, up to 550 GPa and 5000 K. The static lattice 
energy is computed by using the LAPW method, with LDA, PBE, and the recently proposed WC 
functional. The electronic thermal free energy is evaluated using the Mermin functional. The 
vibrational part is computed within the quasi-harmonic approximation using density functional 
perturbation theory and pseudopotentials. Special attention is paid to the influence of the electronic 
temperature to the phonon frequencies. We find that in overall LDA results agree best with the 
experiments. Based on the DFT calculations and the established experimental data, we develop a 
consistent thermal EOS of platinum as a reference for pressure calibration. 



PACS numbers: 64.30.Ef, 05.70.Ce, 71.15.Mb, 63.20.D- 



I. INTRODUCTION 



Platinum (Pt) is a widely used high-pressure standard. Its equation of state (EOS) at 
room temperature has been established by reducing shock Hugonioti'^i'^i^ and by ab initio 
linear-muffin-tin-orbital (LMTO) calculations^ up to 660 GPa. Mao et al^ used the EOS 
developed in Ref. y (Holmes et al.) to calibrate pressure in their compression experiment on 
Fe and Fe-Ni alloy. The bulk moduli measured at the earth's core pressure are substantially 
higher than those extrapolated from seismological observations.-!^ A large pressure offset is 
needed to remove the discrepancy: about 8% at 100 GPa, 15% at 300 GPa. The origin 
of this offset is under investigation. One possibility is the EOS of Holmes et al. seriously 
overestimates pressure.- Singh raised other possibilities.- He noticed that the one-parameter 
EOS of platinum agrees with the EOS of Holmes et al. to 1 % at high pressures, and 
concluded that large systematic error in pressure scale is unlikely. He further proposed that 
the discrepancy is due to the pressure on the sample is different from the one on the pressure 
standard, in the high-pressure X-ray diffraction measurements. 

There are conflicting reports on the uncertainties of Holmes et a/.'s EOS. Dewaele et 
alM- measured the EOS of six metals at ambient temperature to 94 GPa using a diamond 
anvil cell (DAG). By cross-checking different pressure scales they found Holmes et a/.'s EOS 
overestimates pressure by ^ 4 GPa near 100 GPa at room temperature. This conclusion is 
con&med by other groups. While more recent calculations based on density functional 
theory (DFT) suggest Holmes et a/.'s EOS underestimates, rather than overestimates, pres- 
sure. Xiang et alM computed the thermal equation of state of platinum using LMTO and 
a mean field potential method. The pressure they obtained is 5 to 6 % higher than that 
of Holmes et al. at high pressures. Menendez-Proupin et al.~ reached similar conclusion 
using pseudopotentials. Both calculations employ the local density approximation (LDA) 
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However it can also 



functional. And the excess pressure is attributed to LDA in Ref. 
be caused by other factors. In Table II of Ref. [ill, the equilibrium volume decreases as the 
temperature increases. The electronic thermal pressure is negative according to this calcu- 
lation, which is contrary to expectations. Ref. [12 uses an ultrasoft Rappe-Rabe-Kaxiras- 
Joannopoulos pseudopotential from the PWSGF website,— which contains only 5d, 6s and 
6p valence states. Its large cutoff radius (2.6 a.u.) may cause error in studying the highly 
compressed structure. 
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Besides the room temperature isotherm, accurate thermal pressure (-Pth) is needed to 
cahbrate pressure in simultaneous high-pressure and high-temperature experiments. Exper- 
iments cannot easily determine Pth over a wide temperature and volume range.— Conse- 
quently Pth is often estimated by assuming it is linear in temperature and independent of 
volume.-i^ Theory can in principle do better. In quasi-harmonic approximation (QHA), DFT 
calculations give Pth at any particular temperature and volume. It is desirable to combine 
the experimental data with DFT calculations, taking the advantages of both, and construct 
a more accurate thermal EOS for pressure calibration. 

In this paper we have three goals: first is to check the accuracy of the theoretical EOS of 
platinum predicted by different exchange correlation functionals. In contrast with previous 
calculations, we find the room temperature isotherm computed with LDA lies below, and 
nearly parallel to the experimental compression data. The Fermi level of platinum lies in 
the d band and gives a very large density of state (DOS) N{Ep). Its vibrational frequencies 
are more sensitive to the electronic temperature than those of many other metals. A Kohn 
anomaly has been observed in platinum at 90 K.^^ It becomes weaker and finally disappears 
when the temperature increases. Thus our second goal is to discuss the electronic temper- 
ature dependence of vibrations (ETDV) and its influence on the thermal properties. Our 
last goal is to provide an accurate thermal EOS for pressure calibration. For this purpose 
we make corrections to the raw DFT results. We correct the room temperature Gibbs free 
energy G{P, 300 K) to ensure that it reproduces the experimental isotherm, then combine 
it with the DFT calculated temperature dependence to get G{P, T). The thermal EOS and 
thermal properties deduced from the corrected Gibbs free energy are in good agreement with 
the known experimental data. 

II. COMPUTATIONAL METHOD 

The EOS of a material is determined by its Helmholtz free energy P(V, T), which consists 
of three parts: 

F{V,T) = U{V)+F,UV,T)+F,UV,T), (1) 

where U{V) is the static energy of the lattice, F^ihiV, T) is the vibrational free energy, 
Feie{V,T) accounts for the thermal excitation of the electrons. U{V) is calculated by us- 
ing the linearized augmented plane-wave (LAPW) method^ and three different exchange- 
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correlation functionals: LDA, Perdew-Burke-Ernzerhof (PBE),— and Wu-Cohen (WC).— 
The 4/, 5s, 5p, 5d, 6s are described as valence states, others are treated as core states. The 
convergence parameter RK^^^^ is 10.0, and the muffin-tin radius R is 2.08 a.u.. A 16 x 16 x 16 
Monkhorst-Pack^° uniform k-grid is used and the integration over the whole Brillouin zone 
is done with the tetrahedron method. -'^ All the calculations using LAPW are performed with 
and without spin-orbit effect. 

In contrast with the static lattice energy U{V), which is sensitive to the relaxation of the 
core states and requires a full-potential treatment, thermal excitations contribute to much 
smaller energy variations and mostly depend on the valence states. We use pseudopotentials 
to compute the thermal effects. An ultra-soft Vanderbilt pseudopotential^^ is generated from 
the reference atomic configuration 5s'^5p^5d^6s^6p^, including non-linear core corrections.— 
There are two projectors in the s channel, 5s and 6s; two in the p channel, 5p and an 
unbound p at 0.2 Ry above the vacuum level; one in the d channel, 5d. The local component 
is set in the / channel at the vacuum level. The cutoff radii for each channel s, p, d and 
local are 1.8, 1.9, 1.9 and 1.8 a.u., respectively. We use the scalar relativistic approximation 
and spin-orbit effect is not included. This pseudopotential reproduces the LAPW electronic 
band structure, both at the most contracted volume and the GPa experimental volume. 
We find pseudopotentials with different exchange-correlation functionals yield very similar 
electronic band structures for platinum, and we use LDA to compute all the thermal effects. 

We consider 20 different volumes, with lattice constants from 7.8 to 6.2 a.u. ('17.58 to 
8.83 A'^ in volume). For each volume Vi, we use LAPW to compute its static energy U{Vi) 
and the LDA pseudopotental to evaluate its thermal free energy Fvib(Vi,T) and Feie(Vi,T). 
-FVib(Vi,T) is treated within QHA with phonon frequencies dependent on electronic temper- 
ature (denoted as eQHA) as 

fT\V.^ T) = IJ: hu^.,{V.^ r.i.) + k^T Y: ln(l - exp(^^^^:^01^)), (2) 

where Uqj{Vi,Teie) denotes the phonon frequency computed at volume Vi and electronic 
temperature Te\e- In thermal equilibrium the system temperature T, the ionic tempera- 
ture Tion, and Teie are equal. We distinguish these three temperatures to emphasis the 
temperature dependence of phonon frequencies come from different sources. Anharmonic 
phonon-phonon interactions cause phonon frequencies to depend on Tion, but they are omit- 
ted in QHA. Electronic thermal excitations disturb the charge distribution in the crystal 
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and cause phonon frequencies depend on Teie- In the normal QHA used for insulators and 
some metals, this effect is also ignored and u^j has no temperature dependence (except 
through V{T)). Platinum has a larger N{Ep) than many other metals and ETDV may 
have noticeable effects on its thermal properties. To quantitatively measure the influence of 
ETDV, we compare the vibrational free energies at volume Vi and temperature Tj (Tj=500, 
1000, . . . 5000 K) computed without/with ETDV. Without ETDV (normal QHA), phonon 
frequencies are computed at Teie=0 K by using Methfessel-Paxton^i (MP) smearing with 
a smearing parameter of 0.01 Ry. The corresponding vibrational free energy is denoted as 
F^^^{Vi,Tj). With ETDV (eQHA) phonon frequencies have to be computed separately for 
each Tj. This is achieved by using the Mermin functionaP^ and Fermi-Dirac (FD) smear- 
ing. The corresponding vibrational free energy is denoted as F^^^^{Vi,Tj). The difference 
between these two, AFETT)YiVi,,Tj)=F^^^^{Vi,Tj) — F^^^{Vi,,Tj), describes the correction 
caused by ETDV. To get AFetdv at arbitrary temperature between 0-5000 K we fit a 4th 
order polynomial from AFETDv(Vi, Tj) 

AFetdv(V^.,T) = fT^{V.,T)-FX''{V.,T) 

= a,{Vi) ■ T + a2{V,) ■ + a,{V^ ■ + a,{Vi) ■ T\ (3) 

The final vibrational free energy is computed as Fvib(Vi, T)=F^^^(Vi, T) + 
AFetdv(^,T) (we omit the subscript 'eQHA' and denote F^^^ as -FVib). 

Phonon frequencies in the above procedure are determined by density functional per- 
turbation theory (DFPT)^ as implemented in the Quantum ESPRESSO^'' package. The 
dynamical matrices are computed on an 8 x 8 x 8 q-mesh (29 q points in the irreducible 
wedge of the Brillouin Zone). Force constant interpolation is used to calculate phonon fre- 
quencies at arbitrary q vectors. The summation in Eq. ([2]) is evaluated on a 32 x 32 x 32 
q-mesh. 

The electronic free energy FeieiVi-, T) is determined by using the Mermin functional^ and 
Fermi-Dirac smearing. Similar to getting Fvib(V^i,r), we first compute Fgie at every 50 K 
from 50 K to 5000 K, then we fit them to a 4th order polynomial 

Feie(l^„T) = h{V,) ■ T + ■ + h{Vi) ■ + b,{V,) ■ T\ (4) 

Terms other than 62(Vi)T^ represent deviations from the lowest-order Sommerfeld expansion 
Feie = —^N{EF,Vi){kBT)'^, where N{EF,Vi) is the electronic density of states at Fermi 
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energy Ep and volume Vi. We find below 1000 K, keeping only the quadratic term does 
not introduce much error. The infiuence of the higher order terms becomes prominent 
at high temperatures. At 2000 K, the error reaches about 15 %. The fitted quadratic 
coefficient b2{Vi) differs from the Sommerfeld value -fN{EF,Vi)kl by 5 % (^^=8.83 A^) 
to 15 % (1^2=15. 63 A'^). It seems the Sommerfeld expansion works better at high pressures, 
where the electronic bands are more dispersive and N{Ep) is smaller. We combine Feie{Vi, T) 
with the static energy U{Vi) from LAPW and the vibrational free energy Fvib(l^,T) from 
the same pseudopotential to get the total free energy F(Vj,T) at volume Vi. There are 
two popular parameterized forms to fft the total free energy F{V,T), 4th order Birch- 
Murnahan!2^(BM) and Vinet.— We ffnd BM and Vinet are comparable in accuracy to fft 
the static and low temperature free energy, but BM yields much lower residual energies than 
Vinet for the high temperature results. Thus we use 4th order BM to get F{V,T). Other 
thermodynamical properties are computed by ffnite difference. 

All the pseudopotential calculations are carried out with the same plane-wave cutoff of 
40 Ry, charge-density cutoff of 480 Ry, and a shifted 16 x 16 x 16 Monkhorst-Pack mesh. To 
determine the convergence uncertainties of our results, we choose one volume {Vi = 15.095 
A) and recompute its phonon frequencies at T^c = K, with a 24 x 24 x 24 mesh and a 
higher plane- wave cutoff (60 Ry). The two sets of phonon frequencies differ by 0.5 % at 
most. The corresponding differ by 0.07 mRy at 2000 K, 0.18 mRy at 5000 K. The 

inffuence of ETDV is much greater than the convergence uncertainties. For some modes 
phonon frequencies change by 10 % or more from Teie=500 K to Teie=2000 K. The free 
energy correction AFetdv is about 1 mRy at 2000 K. 

III. SUMMARY OF PREVIOUS WORKS 

Besides the two calculationsiii^ mentioned in the introduction, which focus on the ther- 
mal EOS of platinum, there are some other papers related to this subject. Cohen et alM 
computed the static EOS of platinum using LAPW and PBE, and treated it as an example 
to discuss the accuracy of different EOS formations. They found Vinet fitted better than 3rd 
order BM. The accuracy of 4th order BM and Vinet were comparable. Tsuchiya et al.^^ com- 
puted the electronic thermal pressure (Peic) of Au and Pt using LMTO and LDA. At 2200 K, 
Peie is 1.01 GPa for Pt, while only 0.06 GPa for Au. This is caused by the different N{Ep) of 
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the two metals. The small ETDV effect (1 % change in phonon frequency from Teie= K to 
3000 K) observed on gold^^ is consistent with this picture. Wang et al— used LAPW and an 
average potential method to determine the thermal contributions. Then they reduced the 
experimental shock Hugoniot and got the room temperature isotherm of Pt. This isotherm 
is very similar to that of Holmes et al, in spite of the fact that in the latter case, thermal 
pressure is estimated semi-empirically. Ref. |3J computed the static EOS of platinum using 



pseudopotentials with/without spin-orbit effects up to 150 GPa. In the following section, 
we compare our calculations with these previous ones whenever appropriate. 

On the experimental side. The reduced isothermal P-V-T EOS from shock wave experi- 
ments are widely used as primary pressure scales. At present they are also the only experi- 
mental sources for P-V-T data at very high pressures. The shock Hugoniot of platinum was 
first obtained by using chemical explosives.- The reduced room temperature isotherm was 
up to 270 GPa. Holmes et al.^ went to higher compression ratio using a two-stage light-gas 
gun. The final shock Hugoniot is a combination of these two sets of data. In spite of the cru- 
cial role of the reduced shock EOS, its accuracy suffers from low precision in measurements, 
and theoretical simplifications made in the reducing process.—-'-^ With the development of 
DAG and third-generation synchrotron light source, cross-checking different pressure scales 
became feasible. More accurate thermal EOS were obtained by using this method.-ii^ 

Recently, Dorogokupets et alr^ constructed a semi-empirical model to describe the ther- 
mal properties of Al, Au, Gu, Pt, Ta and W, The model contains about 20 parameters, 
which are fitted to the available experimental data on the heat capacity, enthalpy, volume, 
thermal expansivity, bulk modulus and shock Hugoniot. Based on this model they rean- 
alyzed the data in Ref. up to 100 GPa. The resulting EOS, which are consistent with 
the measured thermal properties, are believed to be more accurate than the original in the 
corresponding pressure range.— A simplified version of the model^ yields similar EOS at 
low pressures. However their high pressure extrapolations differ by 2.5 % near 240 GPa. It 
will be interesting to use DFT to explore the EOS at very high pressures, which are still out 
of reach for the current DAG experiments. 
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IV. RESULTS AND DISCUSSIONS 



A. Static Equation of State 

Before studying the EOS at finite temperature, we examine the static EOS computed 
by using different exchange-correlation functionals, and compare them with previous calcu- 
lations. Excluding the thermal effects (which amount to ~ 2 GPa at room temperature) 
helps to identify the origin of their differences. Fig. [1] shows the static pressure vs. volume 
relations using different exchange-correlation functionals. The corresponding EOS param- 
eters are listed in Table [H The experimental data at room temperature are also included 
in the figure to give a rough estimate of the difference. Comparing to the experiments, in 
the entire volume range LDA underestimates pressure while PBE overestimates. WC im- 
proves on PBE, but still overestimates. A detailed comparison between the calculated room 
temperature isotherms (including the thermal effects) and the experimental data will be 
given in Sec. IIVCI DFT has many different implementations, such as LAPW, LMTO, and 
various pseudopotentials. If the calculations are good, they should yield similar results. We 
compare our LDA calculations with previous ones in Fig. [21 Two of our own pseudopotential 
calculations are included for comparison. One is the Vanderbilt pseudopotential that we use 
to compute the thermal effects, denoted as pseudo-1. The other is a Rappe-Rabe-Kaxiras- 
Joannopoulos pseudopotential from the PWSCF website (Pt.pz-nd-rrkjus.UPF), denoted as 
pseudo-2. The static EOS predicted by pseudo-1 is similar to that of LAPW. Their EOS 
parameters differ by no more than 0.5 %. The previous overestimations of pressure are pro 
ably caused by the large cutoff radius and insufficient number of valence e^ 
or another issue related to the negative electronic thermal pressure (Ref. Illl ). 



ectrons (Ref. Il2l). 



Platinum is a heavy element, and its electronic band structure is sensitive to spin-orbit 
effect.— We find inclusion of the spin-orbit effect increases the equilibrium volume, no matter 
which exchange correlation functional is used. This tendency has also been observed by 
Bercegeay et al?'^ in their pseudopotential calculations. However, the EOS parameters are 
not independent of each other. The variation of the equilibrium volume largely compensates 
that of the bulk modulus and the actual pressure difference is within 0.7 % at high pressures. 

Using pseudopotentials instead of the all electron LAPW may introduce error in com- 
puting phonon frequencies, especially at high pressures. Since lattice vibrations are closely 
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TABLE I: Static EOS parameters obtained from LAPW calculations and compared with those 
in literature. Parameters from the pseudopotential calculations (pseudo-1 and 2) are also listed. 
For convenience both Vinet and 4th order BM parameters are shown. Vq denotes the equilibrium 
volume, Kq, Kq and Kq are the isothermal bulk modulus, the first and second derivative of the 



bulk modu 
GPa (Ref. 



us at Vb, respectively. Note their different fitting ranges: 0-550 GPa^ 
, 0-660 GPa (Ref. 
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12), 0-150 GPa (Ref. 



3J), 0-350 GPa (Ref. 



(this study), 0-1000 



3d). Ref. 



34 uses 3rd 



order BM EOS so the corresponding K'J. are not listed. 
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FIG. 1: Static EOS computed by the LAPW method using various exchange correlation functionals. 
(a) Pressure vs. volume curves. Spin-orbit effect is too small to be identified on this scale, (b) 
Pressure difference with respect to the static LDA EOS. In spite of the relative large change 
in EOS parameters, the actual pressure difference with/without spin-orbit effect is very small. 
Experimental data labeled as 'McQueen' are from Ref. [jj, 'Dewaele' from Ref. la, 'Holmes' from 



Ref. 



4. 



related to the force/stress on the atoms, we estimate the error in phonon frequencies by 
analysing the error in static pressure. At high pressures (150-550 GPa), the pressure dif- 
ference between LAPW (with LDA functional) and pseudo-1 is about 1.4 %. The error in 
phonon frequencies caused by using pseudo-1 is likely to be of the same order. Since the in- 
fluence of spin-orbit effect is half of the pseudopotential uncertainty, it is ignored completely 
in the following calculations. 



B. Phonon Dispersion and Its Electronic Temperature Dependence 

Fig. [3] shows the phonon dispersions at the experimental ambient condition lattice con- 
stant a=7.4136 a.u..^ One is computed at Teie=0 K. The other at Toie=2000 K, close 
to platinum's melting point at ambient pressure 2041.3 K.~ The Kohn anomaly (near 
q=[0, 0.35, 0.35]) disappears when the electronic temperature is high, and the vibrational 
DOS varies noticeably. The corresponding corrections to the vibrational free energy, 
AFetdv(^)^)) are shown in Fig. |l](a). AFetdv is always positive. As volume decreases, it 
diminishes and finally becomes negligible. ETDV originates in the thermal excitations of the 
electrons near the Fermi surface, and the number of thermal excited electrons is proportional 
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FIG. 3: (a) Phonon dispersion and (b) vibrational DOS at a=7.4136 a.u.. Solid line corresponds 



to reie=0 K, dashed line corresponds to Teio=2000 K. They are calcu 
pseudopotential. Experimental data labeled as 'Button' are from Ref. \lG 



ated by using the LDA 



to N{Ep) in the lowest-order Sommerfeld expansion. For smaller volumes, the electronic 
bands are more dispersive and N{Ep) decreases, as shown in Fig. |4](b). ETDV diminishes 
accordingly. 

Figure [5] and [6] show the volume thermal expansion coefficient a, heat capacity at con- 
stant pressure Cp, entropy S, and the temperature-dependent part of the Gibbs free energy 
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FIG. 4: (a) Corrections to the vibrational free-energy at various lattice constants, (b) Volume 
dependence of the electronic density of states. 

AG{P,T)=G{P,T)-G{P,T = 300K). Including ETDV removes about half of the discrep- 
ancies between experiments and calculations based on normal QHA. The remaining small 
differences between theory and experiment are attributed to anharmonic phonon-phonon 
interactions^ and electron-phonon interactions."^^ These two effects are of the same order 
of magnitude^ as Fgie for metals, but explicit perturbative calculations to determine their 
magnitudes are computational demanding and beyond the scope of the current paper. We 
notice DFT calculations based on QHA describe well the thermal properties of other metals, 
such as gold,'^^ silver, '^^ copper,— up to melting point. This is in contrast with ionic crystals 
like MgO, where there are large deviations from QHA at high temperatures. It is possible 
that the effects of anharmonic phonon-phonon interaction and electron-phonon interaction 
tend to cancel each other in these metals. Further work is needed to clarify this issue. 

C. Room Temperature Isotherms 

By fitting the total Helmholtz free energy at 300 K, we get the theoretical 300 K isotherms, 
as shown in Fig. [71 Their parameters are listed in Table [Tll In the low pressure range, 
the LDA isotherm and the experimental data are almost parallel. As pressure increases, 
they start to merge. It seems LDA works better at high pressures. Regarding to EOS 
parameters, LDA gives equilibrium volumes closest to the experiments, WC yields closest 
bulk modulus (Kq) and the derivative of the bulk modulus (Kq). Some people^*^ prefer to 
compare pressures from two EOS (labeled as EOS-I and EOS-II) at the same compression. 
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FIG. 5: Thermal expansivity as a function of (a) temperature and (b) pressure. Curves with 
label '(Teie=0)' represent properties computed without ETDV, i.e. computed from F^^^. Curves 
without this specification are the default ones computed with ETDV. '0 GPa corrected' denotes 
the results obtained by adding a phenomenological correction to account for anharmonic phonon- 
phonon interactions and electron-phonon interactions. These corrected thermal data are used to 
construct the final thermal EOS in Table Hill and ITVl The experimental data are from Ref. |46. 



i.e. the same value of V/Vq. Vq is the corresponding equilibrium volume, Vb,i for EOS-I, 
Vo,ii for EOS-II. Such comparisons can give favorable agreement when Kq and Kq of EOS-I 
are close to those of EOS-II, even Vo,i and Vo,ii are quite different.^ As mentioned before, 
the EOS parameters are not independent of each other. It can be fortuitous that Kq and 
Kq agree well. Judged from pressure vs. volume relation, LDA is the optimal functional 
for platinum. It is worth noting the LDA (HGH) pressure vs. volume relation reported 



in Ref. |3j is similar to those obtained in this study. However, Ref. |3j presents data in 
volume vs. compression, and concludes LDA overestimates pressure by 8 GPa near 100 
GPa. In fact, although -R'o.lda (291 GPa from this study) is much larger than -f^cexpo (273.6 
GPa from Ref. |8|), the bulk modulus computed at the experimental equilibrium volume 
Vo.expe (15.095 A^) is 270 GPa, quite close to -ft'o.expe- Thus when plotted in pressure vs. 
volume, the isotherm computed by LDA is nearly parallel with the experimental data in the 
low pressure range. 
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FIG. 6: Thermal properties of platinum, (a) Heat Capacity at constant pressure, (b) Entropy, (c) 
Temperature dependence of the Gibbs free energy at constant pressure. The meanings of the labels 
are the same as above. On this scale it is difficult to discern the improvements on entropy and Gibbs 
free energy, caused by including the phenomenological correction on anharmonic phonon-phonon 
interactions and electron-phonon interactions. Thus those data are not shown. The experimental 
data are from Ref. 47. 



D. Thermal EOS of Platinum for Pressure Calibration 



In the previous sections the thermal properties of platinum is discussed from a pure 
theoretical point of view. We have computed the static lattice energy U (V) using LAPW, 
and found spin-orbit interactions are not important in determining the EOS of platinum. 
We have used QHA to calculate the vibrational free energy F^\\J(y,T\ and found including 
ETDV improves the agreement on the thermal properties. We have calculated the electronic 
free energy Feie(V^, T) using Mermin functional. The resulting thermal properties, e.g. the 
temperature-dependent part of the Gibbs energy AG{P,T), are close to the experimental 
data at GPa, The room temperature isotherm computed by LDA merges to the reduced 
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FIG. 7: 300 K isotherms, (a) pressure vs. volume, (b) pressure difference. 'Dorogokupets' denotes 



EOS from Ref. 



36l . which is very close to the reduced shock EOS (McQueen) near 100 GPa. PBE 



is not plotted as its EOS is way off. 



TABLE II: EOS parameters of the theoretical 300 K isotherms, compared with the experiments. 



Vb,expe is 15.095 A3. Pressure range: 0-550 GPa (this study), 0-660 GPa (Ref. \i and 
GPa (Ref. 5 andQ), 0-270 GPa (Ref. 0). 



12), 0-94 







Vinet 






B-M 








Vo (A3) 


Ko (GPa) 




Vo (A3) 


Ko (GPa) 




(GPa-i) 


LDA 


14.884 


291.25 


5.547 


14.886 


291.65 


5.496 


-0.03232 


LDA (^Cexpe) 




269.96 


5.640 




269.91 


5.626 


-0.03730 


LDA-^ 


15.188 


281 


5.61 










PBE 


15.864 


225.55 


5.751 


15.866 


225.34 


5.741 


-0.04709 


WC 


15.322 


263.93 


5.601 


15.325 


264.72 


5.530 


-0.03580 


Holmes'* 


15.10 


266 


5.81 










Dewaele*^ 


15.095 


273.6 


5.23 










Dorogokupets*^ 


15.095 


276.07 


5.30 










McQueen^ 








15.123 


277.715 


4.821 


-0.01379 



''Reference tl£ 
^Reference |3 ■ 



Reference When Kq is set to 277 GPa, the value measured by ultrasonic experiments, Kq equals 5.08 



GPa. 



'^Reference 



Reference 36|, improved analysis using data from Ref. 



Fitted from the tabulated shock reduced isotherm at 293 K. 
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shock data at high pressures, indicating LDA works better at high pressures. 

Based on these DFT resuhs and all the available experimental data, we try to construct 
a consistent P-V-T EOS of platinum up to 550 GPa and 5000 K. To reach this goal, first 
we need to include the physical effects which are missing in our original model. A phe- 
nomenological term=S AFcorr(y,T) = — |A;b'3'(V/Vo)™T^ is added to the total Helmholtz 
energy to account for the anharmonic phonon-phonon interactions and electron-phonon in- 
teractions, where Vq is the volume of a primitive cell at ambient condition (Vq = 15.095 A). 
The quadratic temperature dependence comes from the lowest order perturbation at high 
temperatures, a and m are two parameters to be fitted. We find setting a equals 10~^ K~^, 
m equals 7 yields good agreement between theory and experiments on a, Cp and AG{P, T) 
at GPa, as illustrated in Figs. [S](a) andlHI^a). The contribution to thermal pressure can be 
estimated by differentiating AFcorr(V^,T) with respect to volume. At 2000 K, APth is 0.38 
GPa when V equals Vq, 0.2 GPa when V equals O.OVq. 

Having obtained accurate AG{P,T), the next step is to get reliable G{P, 300 K). We 
choose the room temperature EOS developed by Dorogokupets et alr^ as our reference 
below 100 GPa. It has been cross checked with other pressure scales, and is likely to be 
more accurate than the reduced shock data of Ref. l| in this pressure range. On the other 
hand, extrapolating an EOS fitted at low pressures to higher range can be dangerous. We 
assume LDA works better at high pressures, and the difference between the exact (obtained 
in an ideal, very accurate experiment) and LDA isotherm approaches to zero as pressure 
increases. 

We compare Dorogokupets's EOS^^ (\/o=15.095 A^, 7^0=276.07 GPa, K'^=h20 in Vinet 
form) with our room temperature isotherm computed by LDA (Vo= 14.884 A^, Ko=291.25 
GPa, /<'q=5.547 in Vinet form). The volume difference between these two at each pressure 
AV{P) = Kxpe,3oo k(-P) — ^LDA,3oo k(-P) IS showu in Fig. [HI Since AV{P) decreases rapidly as 
pressure increases, we use exponentially decaying functions to fit and extrapolate. We correct 
the calculated room temperature Gibbs energy G'lda(-P, 300 K) by setting Gcott{P, 300 K) = 
G'lda(P,300 K) + J^AV{P)dP. The isotherm derived from Gcorr(P,300 K) coincides with 
Dorogokupets' EOS below 100 GPa. the upper limit of their fitting. Above 250 GPa, AV{P) 
is almost zero, and the isotherm derived from Gcott{P, 300 K) is the same as the uncorrected 
one. The uncertainty due to volume extrapolation in the intermediate region (100 GPa to 
250 GPa) is estimated from bulk modulus to be less than 1 GPa. It is worth noting that the 
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FIG. 8: Volume Correction to the theoretical isotherm at 300 K. We use AV{P) = 0.1215 • 
exp(-P/34.0846) + 0.0885 • exp(-(P/109.989)2) to fit the volume difference in the range P < 100 
GPa. The exponential functional form guarantees that it approaches zero at high pressures. It 
happens that this extrapolation agrees well with McQueen's reduced shock data. 



ow 100 GPa, as shown 
will only change the 



established EOS of platinum are quite quite similar to each other be^ 
in Fig. [7](b). Choosing a different reference such as the one in Ref. 
results near 100 GPa by 1 GPa. 

We combine Gcorr{P, 300 K) with the temperature-dependent part of the Gibbs en- 
ergy AG{P,T), and get the corrected Gibbs energy Gcott{P,T) at temperature T. From 
Gcott{P,T) we derive all the other thermodynamical properties. Thermal properties like a, 
Gp, S, which depend on the temperature derivatives of the Gibbs energy, are not affected 
by changing G(P, 300 K). In contrast, the isothermal bulk modulus Kt and adiabatic bulk 
modulus Ks will be influenced, as shown in Fig. [9l 

After corrections, both thermal expansivity and bulk modulus agree with the experiments 
well. We expect the product aKx to be accurate. Integrating aKT we get the thermal 
pressure, p^^(V, T) = P{V, T) - P{V, To) = J^^ aKrdT. The calculated oKt and Pth, before 
and after corrections, are shown in Fig. [TOl Pth(y,T) is often assumed to be independent 
of volume and linear in temperature, i.e. oKt is a constant. Ref. Q assumes the thermal 
energy E(T) = SksT, the thermal Griineisen parameter 7 = 70^/Vo, where 70=2. 4, and 
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FIG. 9: (a) Bulk moduli at GPa. 'Uncorrected' denotes the raw DFT results without any 



correction. As noted in Ref. 



3d . the experimental data in Ref. (4^ is inconsistent. The data points 



we show here are digitized from the graphs in Ref. 



48l (denoted as 'Collard') and 36 (denoted as 



'Dorogokupets') respectively. Ks deduced from the corrected Gibbs free energy agrees well with 



the one computed from empirical models in Ref. 



36l . (b) Bulk moduli as a function of pressure. 



K)=15.123 A3. The thermal pressure is obtained from Mie-Griineisen relation 

E{T)^{V) 3kB^o 



r = 6.57x 10"^ -T (GPa). 



(5) 



V Vo 

In Ref. y, aKx is estimated to be 6.94 x 10^^ GPa/K. Both values lie within the variation 
of the calculated aKx, as shown in Fig. [TOlfa). We find that oKt {Pth) has noticeable 
volume dependence. At fixed temperature, it first decreases, reaches a minimum at about 
y/yo = 0.8, then increases. Such behavior originates in the pressure dependence of the 
thermal expansivity (Fig. [5]^b)) and bulk moduli (Fig. [9]^b)). This feature has also been 
observed in Ref. as shown in Fig. [TO](b) . However it is missing in the previous ab initio 
calculation.— 

Thermal Griineisen parameter 7 = is an important quantity. Empirically it is 

often assumed to be independent of temperature. Its volume dependence is described by a 
parameter q= ^ , and 7 can be represented in g as 7 = 70 . From Mie-Griineisen 

relation, it is obvious that q is related to the volume dependence of oKt- If q equals 1, aK^ is 
independent of volume. If q is greater than 1, oKt gets smaller as volume decreases. In Ref.Q 
q is assumed to equal 1. Fei et al.^ determined 7 by fitting the measured P-V-T data to the 
Mie-Griineisen relation up to 27 GPa. They gave 70=2. 72 and g=0.5. Zha et alr^ extended 
measurements to 80 GPa and 1900 K. Their fit gave 70=2. 75 and g=0.25. Our calculation 
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FIG. 10: (a) temperature derivative of Pthi OiKT=^^^ before/after corrections, (b) Thermal 
pressure Pth at different V/Vq, where Vq is the experimental volume at ambient condition (15.095 
A^). Points are the thermal pressures from Ref. 
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FIG. 11: Thermal Griineisen parameter (a) volume dependence at fixed temperature; (b) the 
corresponding q. 



indicates that the temperature dependence of 7 is small. The volume dependence of 7 is 
shown in Fig. [HI The uncorrected DFT calculation tends to overestimate 7. At ambient 
condition 70 equals 2.87. After corrections, 70=2. 70. The corresponding q equals 2.35, much 
larger than the value obtained in Ref. y and y. We notice previous DFT calculation on 
gold^ also gives much larger q than the value in Ref. y. This is probably due to the small 
pressure range explored in Ref. Isl, and limited number of data points measured in Ref. 10. 



Adding the thermal pressure to the room temperature isotherm, we get the thermal EOS 
of platinum, as shown in Table IIIII We compare our results with two DAG measurements 
in Fig. [T2I Within the error of the experiments, the agreement is reasonably good. For 
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TABLE III: Pressure (in GPa) as a function of compression (I-V/Vq, Vq is the experimental volume 
at ambient condition. 15.095 A'^) and temperature (in K), deduced from the corrected Gibbs free 
energy. Pressures above melting point is not shown. The melting point is determined using 



Tm{P) = 2057 + 27.2 x P - 0.1497 x k up to 70 GPa from Ref. 
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FIG. 12: High temperature isotherms after corrections. Lines correspond to the calculated 
isotherms. The experimental data in (a) are taken from Ref. llSl . except those at P=0 GPa, 



which are obtained by integrating the thermal expansivity listed in Ref. 
Ref. 10. 



461 . Data in (b) are from 



convenience of interpolation, parametric forms of the thermal EOS are listed in Table IIVI 
The P-V-T thermal EOS we obtained are very similar to the one in Ref. 
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below 100 



GPa, This is expectable as we used the 300 K isotherm in Ref. 
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as the reference to 
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TABLE IV: Parametric form of the thermal EOS, P{V)=Po + |Ko[(y/Vb)~^ - (^/Vb)"^] • 
{l + UK'o - mV/Vorl - 1] + IlKoKi; + {K^ - 3){K^ - 4) + f][{V/Vor's - 1]'}. At high 
temperatures, the equihbrium volume may exceed the largest volume we calculate. For better 
accuracy we fit the P-V-T data in three different pressure-temperature intervals: (1) 0-100 GPa 
and 0-2000 K, (2) 50-250 GPa and 0-3000 K, (3) 150-550 GPa and 0-5000 K. Pq denotes the 
starting pressure of the corresponding interval. Vq, Kq^ Kq, and Kq are temperature dependent 
parameters, and are fitted to a 4th order polynomial oq + ait -|- a2t'^ + a^t^ + a^f^, where t=T/1000. 
They have a formal correspondence to the usual 4th order BM EOS parameters, which are defined 
at Po = GPa. 
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correct the room temperature Gibbs energy, and the thermal properties calculated from both 
approaches agree well with the experiments. In this P-T range, the uncertainty of our EOS 
is comparable to the one in Ref. l36| . Above 100 GPa, the uncertainty is about 1.4 %, which 
is difference between the LAPW (LDA) and pseudo-1 static EOS. Other sources of error, e.g. 
convergence uncertainty (0.5 %) and ignoring spin-orbit effect (0.7 %) are smaller effects. To 
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FIG. 13: (a) Theoretical shock Hugoniot compared with the expermental data from McQueen et 
al. (Ref. Ill), (b) Temperature along the Hugoniot. 

check the accuracy of our thermal EOS at high pressures, we start from the corrected Gibbs 
energy and compute the theoretical Hugoniot by solving the Rankine- Hugoniot equation: 

Eh{V,T) ~ E,{Vo,T,) = (P^(V,T) + PoiVo,^^^^, (6) 

where EniVjT), Ph{V,T) are the internal energy, pressure at volume V and temperature 
T. Ei(yo,Ti), PQ{Vo,Ti) are the internal energy, pressure at the initial volume Vq and tem- 
perature Tj. The results are shown in Fig. [131 The predicted Hugoniot pressure is in good 
agreement with measurements. The temperature predicted by DFT is lower than the em- 
pirically deduced value in Ref. ll|. The reduction in Ref. [l neglects the electronic thermal 
pressure, and this may cause overestimating Hugoniot temperature.- 

We end this section by comparing our room temperature isotherm with that of Holmes et 
al.. Below 70 GPa, they are almost identical. At high pressures (200-550 GPa), the pressure 
from our EOS is about 3 % lower than the one from Holmes et al.. Holmes et al. used 
LMTO with the atomic-sphere approximation to get the static EOS. In principle, the full 
potential LAPW method used in this study is more accurate. It seems the EOS of Holmes 
et al. overestimates pressure systematically at high compression ratio. But the magnitude is 
much smaller than the pressure offset needed to compensate the discrepancy between Mao 
et al.^s experiment and seismological extrapolation. The real cause of the discrepancy might 
be a combination of several factors. 
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V. CONCLUSIONS 



In this paper, we report our calculations on the static and thermal EOS of platinum 
using DFT with different exchange correlation functionals. Contrary to previous reports, 
we find the room temperature isotherm computed with LDA lies below, and nearly parallel 
to the experimental compression data. We study the lattice dynamics of platinum within 
QHA, and find the electronic temperature dependence of vibrations plays a noticeable role 
in determining the thermal properties of platinum. Combining the experimental data with 
DFT calculations, we propose a consistent thermal EOS of platinum, up to 550 GPa and 
5000 K, which can be used as a reference for pressure calibration. 
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